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ABSTRACT 

We report a measurement of the streaming motion of the stars in the Galactic bar with 
the Red Clump Giants (RCGs) using the data of the Optical Gravitational Lensing 
Experiment II (OGLE-II). We measure the proper motion of 46,961 stars and divide 
RCGs into bright and faint sub-samples which on average will be closer to the near 
and far side of the bar, respectively. We find that the far-side RCGs (4,979 stars) have 
a proper motion of A < /i >^ 1.5±0.11 masyr -1 toward the negative I relative to the 
near-side RCGs (3,610 stars). This result can be explained by stars in the bar rotating 
around the Galactic center in the same direction as the Sun with Vb ~ 100 kms -1 . In 
the Disc Star (DS) and Red Giant (RG) samples, we do not find significant difference 
between bright and faint sub-samples. For those samples A < /i >^ 0.3±0.14 masyr" 1 
and ~ 0.03 ± 0.14 masyr -1 , respectively. It is likely that the average proper motion 
of RG stars is the same as that of the Galactic center. The proper motion of DSs with 
respect to RGs is ~ 3.3 masyr -1 toward positive I. This value is consistent with the 
expectations for a flat rotation curve and Solar motion with respect to local standard 
of rest. RGs have proper motion approzimately equal to the average of bright and faint 
RCGs, which implies that they are on average near the center of the bar. This pilot 
project demonstrates that OGLE-II data may be used to study streaming motions of 
stars in the Galactic bar. We intend to extend this work to all 49 OGLE-II fields in 
the Galactic bulge region. 

Key words: Galaxy:bulge - Galaxyicenter - Galaxy:kinematics and dynamics- 
G alaxy : s t ruct ur e 



1 INTRODUCTION 

Several groups have carried out gravitational microlensing 
observations toward dense stellar fields, such as the Magel- 
lanic clouds, the Galactic center and disc. Until now, hun- 



distribution is related to the mass function and kinematics 
of the lensing objects. Observed high o ptical depth may be 
explained by the presence of the bar ( KJdatski et al 



Alcock et al 



1997, 2000a gumi et al. 2002). There is sub 
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stantial evidence that the Galaxy has a b ar at its cent er 
Vaucoulcur"s|ll964|; |Blitz fc Spcrgci||l99l|; fctanck et al.| [l994 



1997j ; |Kiraga fc Paezyhski| |1994| ; |Hafner et al.| |200C| ). How- 
ever, the parameters of the bar, e.g., its mass, size, and the 
motion of stars within it, still remain poorly constrained. 
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2001, 2002), and thousands are expected in the upcoming 
years by MOA [], OGLE-III | and other collaborations. 

It is well known that the gravitational microlensing sur- 
vey data is well suited for numero us other scientific projects 
(see Paczyhski 1996; Gould 1996). The studies of the Galac- 
tic structure certainly benefit from this type of data. The 
microlensing optical depth probes the mass density of com- 
pact objects along the line of sight and the event time-scale 
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Stanek et al. ( |1997| ) used the Red Clump Giants (RCGs) 
to constrain the axial ratios and orientation of the Galac- 
tic bar. These stars are the equivalent of the horizontal 
branch stars for a metal-rich population, i.e., relatively low- 
mass core helium burning stars. RCGs in the Galactic bulge 
occupy a distinct re gion in the colour magnitude diagram 
( [gtanek et al, 200C and references therein). The intrinsic 
width of the luminosity distrib ution of RCGs in the Galactic 



bulge is small, about 0.2 mag (Stanek et al 



1997 



Paczyhski 



fc Stanek 1998| ) . Their observed peak and width of the mini 
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nosity function are related to the distance and radial depth 
of the bar. 



Furthermore, Mao & Paczyhski (2002) suggested that 



the proper motion measurements of RCGs in the Galac- 
tic center are useful in constraining the Galactic bar pa- 
rameters. By considering a sub-sample of bright RCGs, one 
should be able to isolate to a sufficient degree the stars that 



are on iverage closer to the near side of the bar. Similarly, 



the stars in a faint sub-sample would be more on the far 
side of the bar. If there is a tangential streaming motion of 
100 kms -1 in the bulge/bar, there should be a detectable 
difference of 1.6 masyr" 1 in the average proper motion be- 
tween the bright and faint RCG sub-samples. Measurements 
of this difference provide constraints on the models of the 
Galactic bar. 

To test the feasibility of the method, in this paper 
we analyze stellar proper motions in one of the fields ob- 
ser ved by the Optical G ravitational Experiment II (OGLE- 



II; Udalski et al. 



|200Cj ). In §| we describe the data. We 
present the analysis method in §^ and results in Discus- 
sion and conclusion are given in §|^. 



2 DATA 

We use the data collected during the second phase of the 
OGLE experiment, between 1997 and 2000. All observations 
were made with the 1.3-m Warsaw telescope located at the 
Las Campanas Observatory, Chile, which is operated by the 
Carnegie Institution of Washington. The "first generation" 
camera has a SITe 2048 x 2048 pixel CCD detector with 
pixel size of 24 /im resulting in 0.417 arcsec/pixel scale. Im- 
ages of the Galactic bulge were taken in drift-scan mode at 
"medium" readout speed with the gain 7.1 e~/ADU and 
readout noise of 6.3 e~ . A single 2048 x 8192 pixel frame 
covers an area of 0.24 x 0.95 deg 2 . Saturation level is about 
55,000 ADU. Details of the instrumen t ation setup can be 
seen in 



Udalski, Kubiak fc Szymahski (1997) 



In this paper we use 266 /-band frames of the BUL_SC1 
field centered at (a,S) 2 oao = (18 ft 02 m 32 s .5, -29°57'41"). 
The time baseline is almost 4 years. There are gaps between 
the observing seasons when the Galactic bulge cannot be 
observed from Las Campanas, each about 3 months long. 
The median seeing is ~ 1.3". We u se the VI photome tric 
maps of standard OGLE template (Udalski et al. 2002) as 
the astrometric and photometric references. 

Only about 70% of the area of the BUL.SC 1 field over- 
laps with the extinction map made by Stanek (1996) and 
is used in the analysis. This ensures that we can accurately 
deredden stellar magnitudes. 



3 ANALYSIS 

The standard OGLE template serves as the fixed astromet- 
ric refe rence in our analysis. In the case of BUL_SC1 field, 



the frame adopted as the OGLE template was taken at JD 



= 2450561.715. In order to treat properly spatial PSF vari- 
ations and frame distortions the field is divided into 256 
subframes before processing. Subframes are 512x128 pixels 
with 14 pixel margin on each side to smooth out transitions 
between the local polynomial fits. The shape of the subframe 



reflects stronger y-axis (declination) gradients due to drift- 
scan mode of observation. The actual proper motion analysis 
includes only 180 subframes, that is 70% of the area of the 
BUL_SC1 field for whi ch the accura te interstellar extinction 
data is available from Stanek ( |1996] ). 

We compute the pixel positions of stars in each of the 
subframes using the DoPHOT package (Schechter, Mateo &i 



Saha 1993). At the start of the processing for each expo- 



sure, the positions of stars in a single subframe are cross- 
referenced with those in the template and the overall frame 
shift is obtained. Using this crude shift we can identify the 
same region of the sky (corresponding to a given subframe 
of the template) throughout the entire sequence of frames. 
For each of the subframes, about ~ 300 brightest (I < 17) 
stars categorized by DoPHOT as isolated are used to derive 
the local transformation between pixel coordinate systems 
of a given exposure and the template. The search radius in 
matching the stars between templete and other frames is 
0.5 pixel. We use first order polynomial to fit the transfor- 
mation. The resulting piece-wise transformation adequately 
converts pixel positions to the reference frame of the tem- 
plate. Typical residuals are at the level of 0.08 pixels. 

An example of time dependence of the position for a star 
with the detectable proper motion is shown in Fig. |lj Also 
presented is the best fit model of proper motion (n a ,ns) 
with and without the differential refraction. The star's co- 
ordinates in the sky at the time t are given as follows: 



a = ceo + + a sin C tan z, 
5 — So + i-ist + a cos C tan z, 



(1) 
(2) 



where a, z and C are the differential refraction coefficient, 
the zenith angle and the angle made by the Zenith, the star 
and the South Pole; oto and So are constants. The parameter 
a is a function of the apparent star colour. Here we neglect 
the parallax effect due to the Earth motion because we are 
interested in stars at the distance of the Galactic center. 

We computed «o, So, fi a , and a for all 46,961 stars 
used to transform coordinate systems (approximately the 
number of used subframes times the typical number of stars 
per subframe, 180 x 300). In cases when the star is measured 
in the overlap region of more than one subframe, the data 
set with the largest number of points is selected. Stars with 
the number of data points fewer than 20 are rejected. 

A sample of fitted proper motions is listed 
in Table hi The complete list of all 46,961 stars 
is available in electric format via anonymous 
ftp from the server astro.princeton.edu, directory 
/sumi/propermotion/buLscl.pm.gz. The list contains 
star ID, number of data points, measured fi a and /is with 
their errors, a, standard deviation (Sdev) of data points 
in the fitting, equatorial coordinates, apparent /-band 
magnitude and V — I colour, and extinction corrected 
magnitude Io and colour (V — I)o estimated using the 
extinction map of Stanek (1996). ID, coordina tes, / and 



V — I for each object are identical with those in Udalski et 



al 



( |2002[ ) 

In Fig. |2j, we present the / and V — I Colour Magnitude 
Diagram (CMD) of stars used in this analysis. We also show 
the correlation between the differential refraction coefficient 
a and the apparent V — / colour (uncorrected for extinction) 
for stars with / < 16. In Fig. pi we plot the uncertainty in fi, 
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Table 1. Sample of fitted parameters. 



ID 


N 


fJl(X U 


% COB 6 
I 1110,0 V I 


^6 

) 




a Sdev 
(mas) 


a 2000 
^ hour ^ 


$2000 
Tripp- "1 




I 


V — I 


Io 
(mag) 


IV — T) 

(V 1) 


11933 


263 


-3.86 


0.49 


5.68 


0.49 


0.70 


9.45 


18.0355861 


-30.3390194 


12 


617 


1.878 


11 


638 


1.241 


11934 


230 


-3.18 


0.52 


4.45 


0.52 


-10.40 


9.37 


18.0359096 


-30.3321479 


12 


254 


0.663 


11 


195 


-0.022 


11935 


263 


1.07 


0.58 


1.95 


0.58 


3.01 


11.34 


18.0348306 


-30.3262084 


12 


747 


2.257 


11 


667 


1.550 


11936 


262 


-4.46 


0.46 


-0.13 


0.45 


34.05 


8.80 


18.0358950 


-30.3248649 


12 


868 


4.805 


11 


790 


4.092 


11937 


243 


-1.24 


0.51 


1.12 


0.51 


31.46 


9.47 


18.0346639 


-30.3174773 


12 


503 


5.206 


11 


400 


4.494 


11938 


246 


1.19 


0.51 


3.11 


0.51 


20.69 


9.53 


18.0349745 


-30.3175164 


12 


492 


4.406 


11 


390 


3.690 


11939 


198 


-4.30 


0.71 


0.76 


0.70 


14.21 


11.82 


18.0351497 


-30.3174968 


12 


247 


3.679 


11 


145 


2.961 


11943 


261 


1.84 


0.54 


3.73 


0.54 


-11.53 


10.50 


18.0349116 


-30.3700535 


13 


220 


0.483 


12 


241 


-0.154 


11945 


262 


-1.43 


0.45 


-3.92 


0.45 


3.06 


8.79 


18.0345371 


-30.3657689 


13 


225 


2.597 


12 


276 


1.979 


11948 


259 


-1.10 


0.41 


0.02 


0.41 


26.13 


7.83 


18.0345576 


-30.3619234 


1.3 


072 


4.608 


12 


129 


3.990 


11950 


266 


0.28 


0.53 


2.26 


0.53 


4.35 


10.33 


18.0367544 


-30.3590652 


1.3 


490 


2.135 


12 


500 


1.509 


11951 


265 


-4.08 


0.53 


1.23 


0.53 


10.62 


10.30 


18.0365689 


-30.3574343 


13 


416 


3.648 


12 


395 


2.999 


11952 


248 


1.59 


0.54 


-6.59 


0.54 


6.22 


10.10 


18.0342311 


-30.3572907 


13 


157 


3.256 


12 


203 


2.626 


11953 


244 


-8.37 


1.42 


-2.66 


1.41 


-3.02 


26.65 


18.0364769 


-30.3556208 


13 


294 


0.671 


12 


243 


0.003 


11955 


265 


-0.56 


0.37 


0.06 


0.37 


4.83 


7.25 


18.0354601 


-30.3520741 


1.3 


700 


2.446 


12 


708 


1.811 



iD=23679, (u ,u ) = (3.B14, -1.396) (mas/yr), a=10.369 (mas/tanz) 
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Figure 1. Time variation of the position in ocos<5 (upper) and 
S (lower) of star ID=23679 (V-I=3.027). Solid line indicates a 
model fit of the proper motion (ii a ,fig) = (3.8, —1.4) (masyr -1 ) 
with a differential refraction a = 10.4 mas(tanz) -1 , and dashed 
line represents the same line without the term of differential re- 
fraction. 



Table 2. Mean of uncertainty in it, < > (masyr 1 ) 
function of Iq (mag), is averaged over Iq i 0.5. 




V-I 



ft 




< °> > 



10 
4.84 



11 

1.40 



12 
0.83 



13 
0.88 



14 
1.18 



15 
1.80 



16 
3.44 



a /j,, as a function of extinction corrected /-band magnitude 
Jo , and we list the mean < (T M > in Table |2| 

To check the measurements we cross-identified our stars 
with the list of high proper motio n obje cts detected pho- 

~ fl200l| ). Out of 74 stars 
are in the region used in 
the present analysis, and 52 were recovered. We plo t the 



to metricallv by Ever fc Wo zniak 
in |Eyer fc Wozniak | ( pOOl]) , 53 



proper moti o ns of those 52 stars as measured in both Eyei 
& Woziiiak (2001 ) (^ew) and this work (zx). The very good 



V-I 

Figure 2. Upper: Colour magnitude diagram of the half of the 
analysed stars. Lower: Correlation between the differential refrac- 
tion coefficient a and the apparent colour (uncorrected for extinc- 
tion) for stars with I < 16. 



correlation gives us certain confidence in our measurements. 
Fig. ^ shows the histogram of /iew (thick solid line) , fi of all 
46,961 sample (dotted line), the sample detected with the 
confidence better than 3<r (dot-dashed line) , better than 5a 
(dashed line) and better than 10cr (thin solid line). We can 
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10 12 14 

I (mag) 

Figure 3. Uncertainty in fj, as a function of the extinction cor- 
rected I-band magnitude Iq. 




Figure 4. Comparison o : 
in both Eyer & Wozniak 



100 



irope r motions for 52 stars measured 
(2001) QitEw) anc l this analysis (fi). 



see that our method is more effectiv e in detecting proper 
motions than frilyer fc Wozniak" (2001). 



4 RESULTS 

For the purpose of proper motion analysis we identify three 
populations of stars: RCGs, Disc Stars (DSs) and Red 
Giants (RGs). Those occupy respectively the upper-right, 
upper-left and lower-right regions in Fig. [| the extinction 
corrected CMD. Within each population we distinguish a 
bright and a faint sub-sample indicated by thick and thin 
boxes in the same figure. Those regions of the CMD are 
defined by the following formula: 



13.5 < Jo < 14.23 (bright RCGs), 
14.39 < I < 14.76 (faint RCGs), 
0.84 <(V- I)o < 1-38 (bothRCGs), 
13.0 < To < 15.0 (bright DSs), 
15.0 < To < 16.0 (faint DSs), 
0.3 < (V-I)o < 0.7 (bothDSs), 
15.4 < T < 15.9 (bright RGs), 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 




100 



jj, (mas/yr) 



Figure 5. Histogram of ^*ew (thick solid line), fj, of all 46,961 
sample (dotted line), that detected with confidence better than 
3<r (dot-dashed line), better than 5<r (dashed line) and better than 
lOcr (thin solid line). 



CO 




(V-Do 



Figure 6. Extinction corrected Colour Magnitude diagram of the 
half of stars used in this analysis. Three populations RCGs, Disc 
Stars (DSs) and Red Giants (RGs) are enclosed by upper-right, 
upper-left and lower-right boxes. Thick and thin boxes represent 
bright and faint samples, respectively. 



16.0 < To < 16.5 (faint RGs), 
0.9 < (V - I) < 1.5 (both RGs). 



(10) 
(11) 



The CMD region for the RCGs was found to maxi- 
mize the significance of the proper motion difference between 
bright and faint samples. 
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Figs. |t] and ^| summarize the main results. In Fig. [7] we 
show a histogram of fi a (left panel) and fis (right panel) 
for bright (thick line) and faint (dashed line) samples of 
all three stellar populations. Fig. shows the contour map 
and the mean of the distribution in \i a and /xs for bright 
(thick line and cross) and faint (thin line and cross) samples 
of RCGs, DSs and RGs. In this figure, the top-left corner is 
toward positive I and bottom-right corner is toward negative 
/ direction. We also show the mean value of the relative 
proper motion < [i > in equatorial and Galactic coordinates 
as well as the corresponding Root Mean Square (RMS) in 
equatorial coordinates for each of the samples in Table ^| We 
assumed that the error in \i of each star is the RMS of the 
related distribution. This RMS is the combination of the 
intrinsic scatter in fj, and the astrometric uncertainty and 
therefore it provides a sensible upper limit for the error in 

< jj, >. In measuring the mean value we rejected high proper 
motion objects with \fj,\ > 20 masyr -1 . Note that the mean 
values < fi > of all sample are not exactly zero, as these 
are not identical to all stars used to align the frames. Stars 
with fewer than 20 data points and detected multiple times 
in the overlapping subframes were rejected. The differences 
are insignificant and we can safely ignore them. 

Finally, we also show the mean positional shift of stars 
in each population as a function of time (JD) in Fig. ^. The 
plots of each samples, i.e., all (dot), bright (filled circle) and 
faint (open circle), are shifted vertically for clarity, because 
only the slope is important. 

In Figs. ^| ^ and Table |§| differences in < fi > be- 
tween bright and faint samples of RCGs are clearly seen. 
The main component is along < m > and reaches A < 
fi >~ 1.5 masyr -1 . The accuracy in A < fi > is about 
a ~ 0.11 masyr -1 , therefore the significance of the differ- 
ence is ~ 14ct. In short, stars at the far side are moving 
toward negative Galactic latitude I relative to the stars at 
the near side. The magnitude and orientation of the proper 
motion of one group relative to the other can be explained 
by the rotation of stars in the bar in the same direction as 
the Solar motion around the Galactic center. The res ult is 
consist ent wi t h A < /x >~ 1. 6 masyr -1 estimated by au 
& Pacz |ynskj J2002[ ) using the [Kiraga fe Paczyhski| ( |l994| )'s 
model with the tangential streaming motion of Vb ~ 100 
kms" . Note that in this analysis the measured < [i > is 
not absolute, but relative. Although the astrometric refer- 
ence frame seems to be close to the one based on stars at 
the near side of the bar, in principle it is not associated with 
any pneumatically defined stellar population. On the other 
hand, the difference A < /i > between the two populations 
of RCGs is properly defined and reliable. 

Both RG samples have same proper motion and values 
of < fj.i > are about a mid between the those of bright 
and faint RCGs. This can be explained by that RGs are on 
average at the center of the bar. 

The disc stars seem to be moving toward positive I with 

< Hi >~ 3.3 masyr -1 with respect to RGs, i.e., the Galactic 
center. This is roughly consistent with < fii >~ 3.2 masyr -1 
estimated assuming that disc stars are on average on the ~ 1 
kpc from the Sun and the flat rotation curve of v r ~ 220 
kms -1 . Here we also assumed that the Sun has an additional 
velocity of 19 kms -1 toward (l,b) = (53°, 25°) relative to 
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the L ocal Standard of Rest (RSL) (Binney 
1987|), i.e., (v Q i,v Qb )= (12 kms -1 , 7 kms -1 ). 



Tremaine 



As for the disc 



Figure 7. Histogram of the distribution of /i a (left) and fig 
(right) for bright (thick line) and faint (dashed line) samples for 
populations RCGs, DSs and RGs. 



stars, their distance may be estimated from their apparent 
brightness. They are typically at the main sequence turn- 
off point for an old population, so very crudely they have 
absolute magnitude Mr,o ~ 4 or so. The bright DSs are at 
about Io = 14 mag, the faint are at about 7o = 15.5 mag. 
So, their distance moduli are approximately 10 mag and 
11.5 mag, which corresponds to the distance from us: 1 kpc 
and 2 kpc, respectively. For the faint DSs, slightly larger 
proper motion with respect to RGs is expected than the 
measured one. This difference might be due to the difference 
in rotation curve along the line of sight. The values < ^ >~ 
—0.3 (bright DSs) or —0.05 (faint DSs) might be explained 
by the solar motion (vQi,VQb) relative to LSR. A detailed 
analysis of these results is beyond the scope of this pilot 
study. 



5 DISCUSSION AND CONCLUSION 

We have measured the proper motion for 46,961 stars in the 
OGLE-II BUL_SC1 field covering Baade's window. We dra- 
matically increased the number of objects with large prope r 
motions in this field compared to Eyer fc Wozniak ( 2001 ) 
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Table 3. Mean values, errors and RMS of fi (masyr 1 ) for star samples discussed in the analysis. 



Population 


N 


< fi a cos8 > 


RMS 


< Ms > 


RMS 


< Mi > 


< Mb > 


all 


46961 


-0.046±0.022 


4.69 


0.091±0.020 


4.39 


0.054±0.021 


0.086±0.021 


bright RCGs 


3610 


0.012±0.060 


3.60 


0.279±0.055 


3.29 


0.245±0.056 


0.135±0.059 


faint RCGs 


4979 


-0.871±0.053 


3.76 


-0.930±0.048 


3.37 


-1.247±0.049 


0.259±0.052 


bright DSs 


2004 


1.669±0.082 


3.67 


2.140±0.084 


3.78 


2.697±0.084 


-0.310±0.083 


faint DSs 


3011 


1.413±0.081 


4.47 


2.229±0.078 


4.26 


2.639±0.079 


-0.045±0.080 


bright RGs 


4911 


-0.501±0.067 


4.68 


-0.356±0.060 


4.24 


-0.565±0.062 


0.242±0.065 


faint RGs 


4770 


-0.516±0.080 


5.55 


-0.378±0.072 


4.97 


-0.592±0.074 


0.244±0.078 



who detected 53 high proper motion objects in the same re- 
gion of the Galactic bulge. We have estimated the difference 
in the proper motion between the bright (near-side) and the 
faint (far-side) samples of RCGs. We found that the far- 
side RCGs have a proper motion of A < n >~ 1.5 ± 0.011 
masyr -1 toward the negative / relative to the near-side 
RCGs. The results fit the picture with stars in the bar rotat- 
ing around the Galactic center in the same direction as the 
Sun. The value A < fx > ^ 1.5 masyr -1 is c onsist ent with 
1.6 masyr -1 estimated by Mao fc Paczynsk] (2002) who as- 
sumed a streaming motion of the bar at vt, ~ 100 fcms -1 . 
The presented method used with the OGLE-II data is sensi- 
tive to the relative streaming motion of stars in the Galactic 
bar down to about A < /i >~ 0.1 masyr -1 . 
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As discussed in the previous section, the measured 
proper motions of DS and RG samples seem consistent with 
the basic understanding of stellar motions in the Galaxy, 
and the Solar motion relative to LSR in particular. This 
consistency supports the reliability of our analysis. 



The question of possible contamination by stars in the 
Galactic disc bears some discussion as we select the samples 
of RCGs using the CMD. Only nearby disc main-sequence 
stars, or evolved disc stars are expected to fall in the same 
CMD region as our RCG samples. The effect can be con- 
s idered neglig i ble b ecause these stars are not very numerous 
(|3tanek et ai] [1994- 



One should keep in mind that all measurements of 
< fj, > presented here are not absolute, but relative to the 
astrometric reference frame which is not well known. This 
problem can be solved by using background quasars that 
would be detected in the near future using the OG LE-II 



variability catalog ( |Wozniak et~al| |2002j ; [Eyer | [2002; ) . The 



primary goal of this paper is to demonstrate that proper 
motions can be measured with very high precision using the 
OGLE-II data, sufficiently high to clearly detect the pres- 
ence of a strong streaming motion (rotation) of stars in the 
Galactic bar. While the reference frame established by the 
large number of stars is not well defined with respect to the 
inertial frame of reference, the relative motions of groups of 
stars are well determined. A detailed analysis of these re- 
sults is beyond the scope of the present study. We intend 
to expand our work to all 49 OGLE-II Galactic bulge fields 
covering a large range of I and b around the bar. A thor- 
ough analysis of all available data is underway and will be 
published elsewhere. 
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Figure 9. The mean positional shift as a function of JD for RCGs 
(Top), DSs (Middle) and RGs (Bottom). The data points for each 
of the samples are shifted vertically for clarity and shown with 
different symbols: all (dot), bright (filled circle) and faint (open 
circle). 
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Figure 8. Contour map and mean of the distribution of fj, a and 
fxg for bright (thick line and cross) and faint (thin line and cross) 
samples within three stellar populations defined in the analysis: 
RCGs (Top), DSs (Middle) and RGs (Bottom). Contours enclose 
30, 60 and 80% of stars. Top-left corner corresponds to positive I 
and bottom-right corner to negative I direction. 
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